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A new mapping between the BFKL equation and Beisert's representation of the XXX Heisen- 
berg ferromagnet with spin s = — | is given. The action of the Hamiltonian operator of a spin 
chain with SL(2) invariance on a symmetric double copy of a harmonic oscillator excited state is 
shown to be identical to the action of the BFKL Hamiltonian on the gluon Green function for the 
azimuthal-angle averaged forward scattering case. A natural mapping between the gluon Green 
function, discretized in virtuality space, and the double harmonic oscillator excited state emerges. 

1 Introduction 

In recent years there has been a growing activity concerning the identification of integrable structures 
in four-dimensional gauge theories. This is mainly due to the interest that this subject has for the 
anti de Sitter / conformal field theory (AdS/CFT) conjecture [ 1 ]. After the seminal works in [2], big 
progress has been made in the mapping of anomalous dimensions of gauge invariant Wilson opera- 
tors in super Yang-Mills (SYM) theory to the spectrum of string theory in different backgrounds. A 
crucial step was to realize that the planar one-loop dilatation operator of TV = 4 SYM maps into the 
Hamiltonian of an integrable spin chain. The problem of finding anomalous dimensions translates 
then into the diagonalization of the corresponding Hamiltonian, and all the techniques developed for 
integrable systems become very useful in accomplishing this task. After those first results, the bet- 
ter understanding of the mapping has allowed more general results for larger orders in perturbation 
theory and for different sectors of the gauge and string theories (for an introduction to the field and 
a wider bibliography see RJ). 

Nonetheless, two-dimensional integrable structures in four-dimensional gauge field theory ap- 
peared well before the AdS/CFT conjecture, in the region of high energy scattering in Quantum 
Chromodynamics (QCD). Non-Abelian gauge theories manifest interesting mathematical proper- 
ties when they are investigated in terms of high energy scattering amplitudes in the Regge limit. 
The most remarkable example [4] is the SL(2,C) invariance present in the impact parameter repre- 
sentation of QCD (and Af = 4 SYM) elastic scattering amplitudes evaluated in multi-Regge kine- 
matics PJ-pl. In this context the Balitsky-Fadin-Kuraev-Lipatov (BFKL) pomeron (with vacuum 
quantum numbers exchanged in the i-channel) can be interpreted as a bound state of two reggeized 
gluons where the Hamiltonian has an interesting operator representation | [T0[ with holomorphic 
separability in coordinate space (JTTJ . The iteration of the BFKL Hamiltonian in the s-channel, 
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describing multiple reggeon exchanges in the generalized leading logarithmic approximation, de- 
fines the Bartels-Kwiecinski-Praszalowicz (BKP) equation fl"3][T4") and was found to have a hidden 
integrability [4 11 15), being equivalent to a periodic spin chain of a XXX Heisenberg ferromag- 
net 1 16-T8). This was the first example of the existence of integrable systems in QCD. A similar 
integrable spin chain, an open one this time, has recently been found in kinematical regions of n- 
point maximally helicity violating (MHV) and planar (N c — > oo) amplitudes in Af = 4 S YM where 
Mandelstam cut contributions are maximally enhanced [ 19]. The importance of Mandelstam cuts in 
the complex angular momentum plane for Af — 4 SYM MHV planar amplitudes was first realized 
in p0l[2T| where corrections to the Bern-Dixon-Smirnov (BDS) iterative ansatz (22) for this class 
of amplitudes were found for the six-point amplitude at two loops. 

There is an interesting connection between the integrable structures appearing in the calculation 
of the anomalous dimension of gauge invariant twist (scaling dimension minus Lorentz spin) two 
operators of spin M in Af = 4 SYM, and in the multi-Regge kinematics. As it was shown in 
(23), the link to the BFKL equation appears upon analytically continuing the anomalous dimension 
function to complex values of M. In particular, the pomeron corresponds to the first singularity 
at M = to — 1, for small cj. The discrepancy between this result and the prediction obtained 
from the asymptotic Bethe ansatz has been subsequently explained in |24| by the calculation of the 
corresponding wrapping corrections for the twist two operators. 

The aim of the present work is to add new hints to the relation between BFKL physics and 
integrable systems. In particular, a new mapping between the two-sites one-loop Hamiltonian of the 
sl(2) sector of the Af = 4 SYM theory, in the double oscillator picture for operators of a given spin 
introduced by Beisert in [25], and a matrix version of the forward BFKL equation, after averaging 
over azimuthal angle, is obtained. 

After this brief Introduction to the subject, in Section 2, a general discussion on the BFKL equa- 
tion is provided, explaining the connection between the non-forward and forward limits. A novel 
discretization in virtuality space is described in detail, highlighting the role and physical interpreta- 
tion of the shift and diagonal operators appearing in this representation. The gluon Green function 
is built in this discrete representation and the limit to the continuum discussed. In Section 3, known 
facts about Beisert's representation of the non-compact SL(2) spin chain are introduced to set the 
ground for a direct comparison with the BFKL equation, in Section 4. The mapping between the 
Hamiltonians in both systems and the gluon Green function in BFKL and diagonal states in the 
SL(2) spin chain is then transparent. Section 5 is devoted to a discussion of the results and to indi- 
cate possible connections with other works. The Conclusions are presented in the final section of 
the paper. 



2 Matrix representation of the forward BFKL equation 

In this work the BFKL Hamiltonian is considered directly in two-dimensional transverse momen- 
tum space (fc), where other components have decoupled into an evolution variable (rapidity Y) 
which plays the role of time (5||9). The four-point amplitude for off-shell reggeized gluons has the 
following momentum flow: 
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With this notation, the BFKL kernel for the pomeron channel has two contributions. The first one 
corresponds to k = 0, i.e., there is no propagator in the s-channel and can be written as 

"Reggeized Propagators" ~ g 2 N c 5^ (qi — q-y J 5^ (q% — q 2 ^ 

d 2 f ^ + f d 2 f ^ ^ . (1) 

r 2 (qi-r) 2 J r 2 (q 2 - r) 2 J 

The second piece has k ^ and corresponds to squaring the Lipatov's vertex, i.e, 

u V ■ ■ » r(2) . - V V\ 9 2 N C (qi^ 2 + ftjVl' 2 . 

■■Emission" ~ ^ j (g x + q 2 - g t - q 2 ) _ 2 ^, 2 p- (91 + q 2 ) . (2) 

After a Fourier transform of these two expressions, and complexifying the transverse momenta, 
Lipatov found the very interesting SL(2,C) invariance of this Hamiltonian [4|. In the present work, 
however, the focus lies on the forward limit, with zero momentum transfer q — 0. It is noteworthy 
that in this case the contribution from the "Reggeized Propagators" reads 

"Reggeized Propagators" ~ 2g 2 N c 5 {2) (q x - qi ) d 2 r — — -, (3) 

V J J r 2 (q 1 -r) 2 

while the "Emission" piece simplifies to 

g 2 N c g 2 N c 

"Emission" ~ e T-^ = 2—^ -j-^. (4) 

k 2 (<7i ~ q 1 ) 2 

It is in this forward case that it truly represents a real emission since now the amplitude corresponds 
to the 2 — > 3 inelastic process. A further simplification is very convenient: to integrate over the 
azimuthal angle formed by the two transverse momenta qx and q x . 

Once this is done the BFKL equation for forward scattering can be cast in the simple form 

d<p(Q 2 ,Y) f°° dq 2 f a 2iriin( g 2 ,Q 2 ) 2 \ 

where a = a s N c /n, the integration takes place over the gluon virtuality q 2 (= q 2 ) and the cor- 
respondence with the previous notation is q 2 = Q 2 and q x 2 = Qq. The term with <p(Q 2 ,Y) 
corresponds to Eq. ^ and the one with (p(q 2 , Y) to Eq. Since the forward limit has been taken, 
(p(Q 2 , Y) is the cut reggeized gluon four-point function for a given rapidity Y with the initial con- 
dition ip(Q 2 , Y — 0) ~ S(Q 2 — Qq), and it corresponds to the sum of the squares of the 2 — ^ 2 -\- 71 
emissions amplitude over any number n of real gluon emissions. 
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To find the gluon Green function it is convenient to write Eq. Q in the form 

d<p{Q 2 ,Y) _ f 1 dx ( 2 1 



<p(xQ\Y) + -p[^ 1 Y) -2<p(Q 2 ,Y)\, (6) 
Oiui Jq ± x I x \ x 



and then introduce a Mellin transform to obtain 



e ^x(7) ; ( 7 ) 



*L (9- 

2m \Ql 



<p(Q 2 ,Y) = f 

J a 

X(7) = 2^»(1) - ^(7) - ^(1 - 7), (8) 



with tp being the digamma function and < a < 1. It is well known that, for asymptotically large 
values of the rapidity variable Y, this integral tends to 



aY4loe2- 
e B 



tp(Q z ,Y) ~ , (9) 

with t = log Q 2 /Qq, This implies the following diffusion equation for the function cf> = ipe 1 ! 2 j%: 



■2, 



- 41 ^ 2 ^ 14 C(3)^, (10) 

which shows that there exist two different flows for the virtualities of the i-channel gluons in the 
BFKL ladder: one towards the infrared (IR) and one towards the ultraviolet (UV). These IR/UV 
flows are symmetric since the eigenvalue function x(7) i s invariant under 7 -H- 1 — 7. 

It is possible to gain further insight into this picture. For clarity, it is convenient to mix Eq. ([5]) 
and Eq. (|6|l using the representation 



dip{Q 2 ,Y) [°° df 

adY .Ln2 \l 2 \ 



<p(l 2 + Q z ,Y) + nv(Q 2 ,Y), (11) 

Q 2 ' 



where 



n = -2 f (12) 

Jo l ~ x 

The integration region with I 2 < corresponds to a virtuality flow into the IR, and that with I 2 > 
into the UV. The term proportional to il is diagonal in momentum space since it does not change the 
virtuality of the original gluon. It is formally infinite and allows for the regularization of the infrared 
divergence in the integration of the other term. This separation is introduced because it will make 
explicit interesting properties of the Hamiltonian (or evolution kernel since this is the forward case). 

The space of virtualities can be discretized using I 2 — n A, Q 2 — N A, dl 2 = A and the 
notation <fr n = tp(l 2 , Y). The discrete version of Eq. 
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then reads 



—^7 = > h> +S2 0JV. (13) 

n—l n—1 

Note that, due to scale invariance, there is no dependence on the "virtuality unit" A. 
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To find a matrix representation of the kernel it is useful to introduce the TV -dimensional vector 
= (0i 7 02 1 • ■ ■ i <?W)*> the extended oo-dimensional vector = (tfii, <j>2, ■ ■ ■ , 4>n, ■ ■ . )' and write 



Eq. ( 13 1 in the form 



d(j) 
adY 



H 



N ■ 9oo, 



(14) 



with the kernel being the following semi-infinite matrix with N rows and oo columns: 



1 



n N = 



i 
n 
i 



V N-l 

In terms of components this is equivalent to 



l 

N-2 



1111 

MM 

2 3 4 

"ill 



1 n 1 



(15) 



iV-l tf+n 

E — 



00 & 



(16) 



n— 1 n— 1 

for 1 < i < JV, 1 < j < oo. The virtuality flow towards the IR and UV can be described using shift 
operators: 



Sir 
(<Suv) 



hi 



= 4 + v 



More formally, the Hamiltonian reads 

JV-l 

Hn — 

n=l 



E (^ir)" , (Suv) T 



n=l 



J2 ^ - log (l - Sir) - log (l - Suv) + O i. 

i=JV 



(17) 
(18) 

(19) 
(20) 



In this framework the diffusion picture is clear. At Y = the evolution starts with the initial 
condition ip(Q 2 , Y = 0) ~ 8(Q 2 — Qq), where Qq = Nq A. This corresponds to a single entry in 
the initial condition vector, i.e. 



Putt,. :,<&,■ 



with 



A ' 



(21) 



(22) 



As each emission takes place at a rapidity strongly ordered with respect to the previous ones, it is 

(23) 



possible to write the solution to Eq. ( 14 1 using the following formal expression 



= + dyi ( Q! ^ JV ) + J dyi^aUN^j J dy 2 (ctH 
+ J dyi {ctH N ^j J dy 2 (aH N ^j J dy 3 (aH 



(24) 
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The action of the shift operators in the kernel desplaces the original single "cell" in Eq. ( |22| i towards 
lower or higher entries in the vector. This corresponds to the emission of real gluons in the s-channel, 
which modify the virtuality of the reggeized gluons "propagating" in the t-channel. The diagonal 
piece in the kernel, f2 1, does not modify the virtuality and corresponds to the "propagation" of 
a reggeized gluon with constant q 2 . O formally being an infinite quantity corresponds to a non- 
infrared-finite gluon Regge trajectory. Note that the solution to the original Eq. (JTTJ is given by the 
last component of <f> in the continuum limit where N — > oo, A — > 0, while keeping NA = Q 2 fixed. 



3 sl(2) spin chain Hamiltonian 

As it was mentioned in the Introduction, there has been big progress in the understanding of the 
planar limit of the M = 4 S YM dilatation operator by mapping it to the Hamiltonian of an integrable 
one-dimensional spin chain. The part of the theory of interest in this work is the non-compact 
bosonic sl(2) closed subsector where states are constructed with scalar fields <£> (SO(6) Yang-Mills 
bosons) and their spacetime covariant derivatives P$, which scale under the sl(2) subgroup of the 
Lorentz group. The corresponding charges are related to the SO(4,2) group and to the SO(6) 1Z- 
symmetry. 

The commutation relations of the sl(2) subalgebra of the superconformal algebra are 
[JW,JH] = -2.7(3), [JO), J(±)] = ±J(±). The spin chain Hamiltonian representing the one- 
loop anomalous dimensions of Af = 4 S YM operators with spin S — 1 in the planar limit of the sl(2) 
sector reads 

Hf^e(S - N)(a\) N - 1 (4) s - N \00) = 

0(S-l)(a\) l - 1 (al) s - l \OO),(25) 



1=1 



where A = is the coupling and h(N) = YliLi J = i>(N + 1) — "0(1) is the harmonic number. 
The notation in terms of harmonic oscillators (a^)™|0) = ^(2?)™$ (with a|0) = 0, [a,a^] = 1), 
which corresponds to a site in a one-dimensional lattice (the total number of lattice sites is equal 
to the total TvL-charge), has been used, n indicates that a given site is in the n-th excited state with 
respect to the Tr(<& 2 ) vacuum. These excitations are classified in the spin s = — ^ representation of 
sl(2), which is infinite-dimensional. There is an overall trace in the operator cyclically ordering the 



different sites in the spin chain. Eq. (25 i was introduced by Beisert in [25] and it corresponds to the 
nearest-neighbor one-loop Hamiltonian of an integrable XXX spin s = — \ chain. More explicitly, 
this Hamiltonian is invariant under the sl(2) generators 

J iV = a i(! + a\ a i) +4(! + 4 a 2), (26) 

Jl2 = «1 + a 2, (27) 

j[f = l + 4a x +4«2- (28) 
It is also interesting its action on the following state: 

oo i 

< 2 2) (4-4) P |oo) = E(- 1 ) p " s w^v <2 2) g(p-^)(4) s (4) p ~1oo) 

s=0 >' 

= 2h(p) (a\ -4) P |00). (29) 
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4 The hidden BFKL / XXX 1 spin chain mapping 

To establish the mapping between the BFKL equation and the Hamiltonian associated to the dilata- 
tion operator for the sl(2) sector of the M = 4 S YM theory it is needed to treat the diagonal operator 
£1 in Eq. ( 12 1 of the BFKL kernel in a different way. The discretization of Eq. (Hi can be also 
written in the form 

n=l V 7 n=N+l v 

JV— 1 , oo , 

= E -%-2h(N-l)<f> N , (30) 

n=l n=N+l 

This is a valid representation up to O (^jf-^J terms, which are negligible at large N, where the 
mapping will be defined. The difference between this representation and the previous one is in the 
regularization of the infrared divergence present in the gluon Regge trajectory. Different regulariza- 
tions can be used which are equivalent in the N — > oo limit. 
The new representation of the BFKL Hamiltonian now reads 

oo A n 

= -E — -log(l-<Sm)-log(l-<Suv)+^ (31) 

n=N 

where 

(g) = -2h(i-l)Sj. (32) 
^ ' i,j 

To connect with the XXX_ i spin chain Hamiltonian acting on a state represented by a double 
copy of a harmonic oscillator it is needed to do this on a diagonal state with the same number of 
derivatives in each oscillator (both are in the same (N — l)-th excited state): 

Hf^eiN-l^atf-Haif-'m = 

- A E (jZ^Y - 2h ( N ~ 1 ^) 0(2N-l-l)(afy- 1 (a+) 2N - 1 - l \00), (33) 



which is obtained by setting S = 2N — 1 in Eq. (25 1. 

The equivalent Hamiltonian for the leading-order BFKL calculation in QCD and TV = 4 SYM 
theory can be found by simply rewriting the right hand side of Eq. < (30] > in the form 

U^<t>N = ^(^^-2h{N-l)8?y i . (34) 



From the comparison of expressions ( 33 1 and ( 34 1 a natural mapping between the discretized BFKL 
gluon Green function and the double harmonic oscillator state emerges: 

4>i O 9(2N -l-l)(a+) l ^ 1 (a+) 2N - 1 - l \00). (35) 
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To study how the convergence to the continuum N — > oo limit takes place it is useful to work with 
the function 



1=1 

(N-l oo \ 

E- E 
Z=l i=JV+l/ 



N — I \ I 



NV ^2 



E 



(36) 



In FigjTJit is numerically shown that limjv-j-oo Xw (7) = x(l) °f Eq. ([I]) in the range < 7 < 1. 
The convergence in N is not uniform in this region since it is much faster for small values of 7. 
Analytically, the continuum N — >• 00 limit can be found using N = 1 /A and / = xj A with A — >• 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



0.9 1 



Figure 1: The function xn{j) coincides with the BFKL eigenvalue x{l) for 7 ^ already with 
a low number of excitations N while it requires a very large number of them to match the 7 ~ 1 
region. 



to obtain 



1 dx 



(7) - / , (x^+^- 1 -2) = X ( 7 ). 

Jo 1 — x 



(37) 



5 Discussion 



In this Section possible consequences and open questions for the mapping between the BFKL gluon 
Green function and the double harmonic oscillator state are discussed in some detail. The interpre- 
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tation and answers to these are the bulk of current investigations. 



5.1 Anomalous dimensions and wrapping contributions 

The main assumptions in the asymptotic Bethe ansatz (ABA) were all-loop integrability and the 
existence of PSU(2, 2|4) symmetry in the S-matrix. The ABA breaks down at order <? 8 (four loops) 
for twist two operators due to wrapping corrections. The mapping proposed in this work should shed 
some light on the calculation of wrapping corrections in the spin chain context. A simple procedure 
would be to take the N — >• oo limit, corresponding to an infinite number of magnons, in the action of 
the sl(2) Hamiltonian on the states in the chain and repeat the steps which in the BFKL approach lead 
to the calculation of the leading O ((a/uj) n ) contributions to the anomalous dimension of twist-two 
operators with spin M = u> — 1, for uj — > 0. To follow this path, Eq. (|7jl is written in the form 



2ni J 2iri \QlJ u - aXooil) ' 

where Xoail) — 2^(1) — i'il) — ^(1 — 7)> an d then notice that it is more convenient to perform 
first the integration in 7, i.e. 

w = a-Xoohcu), (40) 

in order to extract the dominant contributions in the collinear region, which corresponds to the 
Q 2 ^S> Qq limit. Eq. (40 » then generates the anomalous dimension: 



(41) 



if the expansion 

w - a I — + 2 £ C(2i + 1) 7 ' L J (42) 

\ 7w L=l / 

of Xoo (7w) is used. It is only in the N — >• 00 limit that the complete expansion in terms of Riemann 
£ functions appears and the correct wrapping terms in the anomalous dimension are found. 

The common understanding is that it is in the long-range interactions of the spin chain where a 
mapping to higher-loop results in the gauge theory side should appear. It is in this sense that it would 
be very interesting to identify the correct variable in the sl(2) spin chain dual to the rapidity Y in the 
BFKL framework and around which to generate a similar exponentiation of the Hamiltonian as in 



Eq. (38 1 and its asymptotic expansion at large Y in Eq. d9). This might lead to a simple method to 
obtain wrapping terms in the spin chain which could be also applied to Hamiltonians associated to 
other symmetries. 

5.2 NLO kernel and comparison with other approaches 

The mapping here discussed might also hold at higher loops both in the sl(2) spin chain and in 



the BFKL formalism. The extension of Eq. (34i to next-to-leading order, i.e., the discretization in 
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virtuality space of the azimuthal-angle averaged forward next-to-leading order BFKL equation in 
Af = 4 SYM can be written in the form 

t/BFKL 
rt NLO 9N 



(43) 




N J 2(N + l) 



It will be worth to investigate the interplay between this formula and any higher order representation 
of the sl(2) spin chain in future works, together with its interpretation in terms of wrapping contri- 
butions (the subleading in uj corrections for the anomalous dimension of Eq. ( pT) are found using 
Eq. (43 1). Note that Eq. (43 i would still correspond to a two-site spin chain with two states in their 
N — 1 excited level. At this point it is fair to bring into the discussion the results found by Shuvaev 
and Wallon in p7| where an expression as in Eq. |5]) was found with (p corresponding there to the 
distribution of soft photons in a charged source, with no relation to high energy scattering in the 
Regge limit or the BFKL equation. 

From a more mathematical point of view, it is likely that a link can be found between the matrix 
representation of the sl(2) Hamiltonian here unveiled and the work of Adler and Van Moerbeke 
in [28 1, where it was investigated how the discrete sinh-Gordon equation leads to the Toeplitz 2- 
Toda lattice. This new lattice has r-functions which are annihilated by operators living in a SL(2,Z) 
subalgebra of the Virasoro algebra, and have a very similar structure to Eq. < fT~9| > if the harmonic 
weights 1/n are identified with time variables in the IR/UV directions. It is on these time variables 
where the Virasoro operators act. 

To conclude, it is noteworthy that in [29] a coherent state representation for the Hamiltonian 
of the spin chain with sl(2) symmetry was derived. The action of the Hamiltonian on the coherent 
states rii , n 2 then reads 



(n 1 ,n 2 \'Hf§' \n ll n 2 ) 



= log 1 



(Hi - n 2 ) 
4 



(44) 



with ^ living on a two-dimensional hyperboloid for sl(2). It will be interesting to study the relation 



with the Eq. ( 3 1 1 derived in this work. 



6 Conclusions 

We have unveiled a new mapping between the BFKL equation and the s = — \ XXX spin chain. This 
mapping requires of a simplification of the original BFKL equation to its forward limit and averaging 
over the azimuthal angle dependence of its kernel. The remaining physical variable corresponds to 
the virtuality of exchanged gluons. Discretization in this variable allows for a direct link between 
the gluon Green function, i.e., the solution to the BFKL equation, and the two harmonic oscillators 
state representing the two-site sl(2) spin chain. This is the main result of our analysis and it is 
presented in Eqs. ( |33|34|35] >. A new matrix representation of the BFKL equation has been given in 



Eq. ( 15 1, together with its physical interpretation in terms of a diffusion process into infrared and 
ultraviolet regions of virtuality space. We have proposed some directions for future research in the 
interpretation and better understanding of this mapping. These include the analysis of higher order 
corrections, the relation to Toeplitz lattices, the study of the mapping in terms of coherent states and 
how it could be related to the calculation of wrapping effects. 
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